Obtaining information on short-range correlations from inclusive electron scattering 
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In view of recent data from the Thomas Jefferson National Accelerator Facility (JLab) on in- 
clusive electron scattering off nuclei at high momentum transfer {Q^ > IGeV^) and their current 
analysis, it is shown that, if the scaling variable is properly chosen, the analysis in terms of scaling 
functions can provide useful information on short-range correlations (SRC). This is demonstrated 
by introducing a new relativistic scaling variable that incorporates the momentum dependence of 
the excitation energy of the {A — 1) system, with the resulting scaling function being closely related 
to the longitudinal momentum distributions. 
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Obtaining information on short-range correlations 
SRC in nuclei is a primary goal of modern nuclear physics 
Interest in SRC stems not only from the necessity 
to firmly establish the limits of validity of the standard 
model of nuclei but also from the impact that the knowl- 
edge of the detailed mechanism of SRC would have in 
understanding the role played by quark degrees of free- 
dom in hadronic matter and the properties of the latter 
in dense configurations 0. Recently, evidence of SRC 
has been provided by new experimental data on inclu- 
sive [A{e, e')X] @, 111 and exclusive [A{e, e'pN)X and 
A{p,pN)X] lepton and hadron scattering off nuclei at 
high momentum transfer {Q^ > 1 GeV"^) (see Ref. Q 
and references therein quoted). In inclusive scattering 
the observation of a scaling behavior of the ratio of the 
cross section on heavy nuclei to that on the deuteron Q , 
for values of the Bjorken scaling variable 1.4 < xbj ^ 2, 
and to that on '^He [1], for 2 ^ xb ^ 3, has been in- 
terpreted as evidence that the electron probes two- and 
three-nucleon correlations in complex nuclei similar to 
the ones occurring in two- and three-nucleon systems Q . 
It should be pointed out, however, that whereas exclusive 
processes can directly access the relative and center-of- 
mass motions of a correlated pair in a nucleus 0] , obtain- 
ing information on these quantities from inclusive scat- 
tering is, in principle, more difficult. Various approaches 
based on scaling concepts have therefore been proposed, 
going from the already mentioned scaling behavior of the 
cross section ratio plotted versus x bj , to the scaling be- 



havior of the ratio of the nuclear to the nucleon cross 
sections plotted versus proper scaling variables; among 
the latter, a process that has been most investigated in 
the past is the so-called Y-scaling, for it is believed that 
this may represent a powerful tool to extract the high- 
momentum part of the nucleon momentum distribution 
which is governed by SRC [!,[§]. It is the aim of this 
Rapid Communication to critically reanalyze the concept 
of Y-scaling, mainly because of i) the lack of a general 
consensus about the usefulness of such a concept and 
ii) a strong renewal of interest in Y-scaling owing to re- 
cent experimental data on A{e,e')X reactions from the 
Thomas Jefferson National Accelerator Facility (JLab) 
13, [H. We will show that the analysis of inclusive scat- 



tering in terms of proper Y-scaling variables could indeed 
provide useful information on SRC; to this end, follow- 
ing the suggestion of Refs. 13, 13, 14 1 a new approach 



to Y-scaling and its usefulness will be illustrated in de- 
tail. Let us consider a virtual photon of high momentum 
impinging on a nucleus A (with mass Ma) and knocking 
out, in a quasielastic process, a nucleon N (with mass 
ttin) having momentum k = |k| and removal energy E. 
The latter is defined as the energy necessary to remove 
the nucleon from A leaving the residual nucleus {A — 1) 
(with mass Ma-i) with intrinsic excitation energy E\_-^ 
(i.e., E^mN + Ma-i -Ma + E\_^ = S„„, + E\^^). 
In Plane Wave Impulse Approximation (PWIA) and us- 
ing the instant form of dynamics, the quasielastic cross- 
section reads as follows 
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with energy conservation {Mj^_^ — Ma-i + E^_-^) 



ly + MA^ Jm% + ik + q)^ + J MXU+'^'' (2) 



and momentum conservation q = p + Pa-i- Here 
1^ — 61 — 62 and q = ki — k2 are the energy and three- 



momentum transfers (Q^ ~ — = 4ei €2 sin^ | with 
q = |q|), aeN is the elastic electron cross section off a 
moving off-shell nucleon with momentum k = |k| and 
removal energy E, and P^{k,E) is the spectral func- 
tion (normalized to one) of nucleon N (i.e., the joint 
probability to have a nucleon with momentum k and 
removal energy E); eventually, p and Pa-i arc the 
momenta of the undetected struck nucleon and the fi- 
nal {A — 1) system. Considering, for ease of presenta- 
tion, isoscalar nuclei, one has P^{k,E) = Pj^{k,E) — 
P,f{k,E) = P^{k,E) = P^{k,E) + P^{k,E), where 

I 



where SeN is the electron-nucleon cross section integrated 
over the polar angle, and F^{q, v) is the nuclear structure 
function 

F^{q,i^) = 2Tr J dE J k dk P'^ik, E) {A) 

Eq. ([3]) is obtained by eliminating the 5 function by in- 
tegrating over cos a = (k • q/kq), with the limits of inte- 
gration resulting from the condition — 1 < cos a < 1. We 
can now introduce a scaling variable Y — Y{q^ v), which 
is only required to be a function of q and v (and any 
arbitrary constant) so that, no matter with the specific 
form of Y, the cross section and the structure function 
can be expressed not in terms of the two canonical inde- 
pendent variables q and v, but, without loss of generality, 
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where /"^(Y) = 27r kdkn^(k) represents the longi- 
tudinal momentum distribution^ and 

00 fc„„„(g,y,£;) 
S'^((7,y) 27r j dE j kdtP^{k,E) (6) 

-E™„. \Y\ 

is the binding correction ,^], which, through 

kminiq,Y, E), is governed by the continuum energy 

spectrum of the final (^ — 1) system, unlike f"^{Y), 



P^ik, E) = {l/A) E„6F Aa,n^{k) 5{E - e<,) is the (triv- 
ial) shell-model part [with Aa denoting the occupation 
number of the single-particle state a with removal en- 
ergy Ca and momentum distribution na{kl)\, and Pi is 
the (interesting) part generated by NN correlations. The 
spectral function is linked to the momentum distributions 
by the momentum sum rule n'^{k) = J P^{k, E) dE = 
J P^ik, E)dE + J P^{k, E) dE = n^{k) + n^{k). It has 
been shown Q that at high values of momentum trans- 
fer, after integrating over the direction of k, Eq. ([T]) can 
be written, to a good approximation, as follows: 



(3) 
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in terms of q and Y — Y(q, i>). Correspondingly, a scal- 
ing function F'^{q,Y) is introduced; this is nothing but 
Eq. (HI) with v replaced everywhere by Y; if, under cer- 
tain conditions, F^{q,Y) F^{Y), Y-scaling is said to 
occur and, depending on the physical meaning of Y and 
F^(Y), various information on nucleons in nuclei could 
be obtained. To simplify our analysis, let us consider high 
values of the momentum transfer, when EmaxiQTY) and 
kmax{<l, Y, E) become so large that, because of the rapid 
falloff fall-off of P^{k, E), they can be replaced by 00 (al- 
though in actual calculations we use the correct values of 
these quantities); in this case, the q and v dependence 
of the scaling function is governed only by fcmm (9, Y^ E) , 
and it is trivial to show that, by adding and subtracting 
a proper term, the scaling function can be cast in the 
following general form 



(5) 



which is integrated over all excited states of {A — 1). 
The quantities f^{Y) and n"^(fc) are linked by the 
relation n^{k) = [df^{Y) / dY]/[2TTY],k = so that 
if f^(Y) could be extracted from the experimental 
data, n^{k) could be determined. Unfortunately, 
such an extraction is hindered by the presence of 
7^ depending upon the difference between Y 
and fc„iin and therefore upon the definition of the 
former. The binding correction is absent only in 
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the deuteron, since E = Emin = cost = 2.22 MeV , 
so that Y kmin{q,v, Emm), {q,Y) = and 
F^{q,Y) = f^{Y). The final state interaction (FSI) of 
the struck nucleon invahdates the PWIA, but in spite of 
that, an approach was developed in the past to reduce 
the effects from both the binding corrections and FSI [8j ; 
the approach is based upon the widely used relativistic 
scaling variable Y = y which is obtained 

by setting in the energy conservation equation [Eq. ^] 
k — y, (\/kq = 1 and, most importantly, E\_^ — 0; 
thus y represents the minimum longitudinal momentum 
of a nucleon having the minimum value of the removal 
energy E = Emin- In the asymptotic limit {q — > oo), 
Eq. ([5]) scales in y and becomes the asymptotic scaling 
function F^{y) = f^iy) - B^{y), that is Eq. (P with 

Y and kmin{q,Y,E) replaced by y and k^^^{y,E), re- 
spectively (scaling in this variable also occurring within 
a relativistic description of the deuteron Unfor- 
tunately, owing to the presence of B^{y), F^{y) is not 
related to a momentum distribution so that, in principle, 
the experimental longitudinal momentum distribution 
fexiy) ciiid, consequently, nf^{k), cannot be extracted 
from the data. Let us briefly recall how this problem 
was addressed in Ref. fsj. The experimental scaling 
function F^^{q,Y) = (jl^^{q,Y)/{[Zsep{q,v,k,E) + 
Nsen{q,v,k,E)]{Eplq)}(^i,^^^^^E^.^) exhibits, when 

Y = y, & strong q dependence owing to the FSI and 
binding effects and differs from the asymptotic scal- 
ing function F^{y\ The latter, however, has been 
obtained in Ref. [8( by extrapolating to g ^ oo the 
available values of F^{q,y), on the basis that FSI can 
be represented as a power series in l/g and dies out at 
large q^, a conclusion that has been reached by various 
authors (see e.g. Refs. The experimental longitu- 
dinal momentum distribution f^{y) has thereby been 
obtained by adding to F^{y) the binding correction 
B^{y) evaluated theoretically, and n^^{k) has been 
obtained by n^(fc) = -d[F^{y) + B ^ (y) / dy] / [2t: y], 
k = \y\. Such a procedure affects the final results 
in terms of large errors on the extracted momentum 
distributions, particularly at large values of fc; in spite 
of these errors, the extracted momentum distributions 
at fc > 1.5 — 2 fm^^ turned out to be larger by orders of 
magnitude from the prediction of mean-field approaches, 
and in qualitative agreement with realistic many-body 
calculations that include SRC. To make the extraction 
of f^{y) as independent as possible from theoretical 
binding corrections, in Ref. ^12] another scaling variable 

Y — yew has been introduced; this scaling variable 
incorporates relevant physical dynamical effects left 
out in the definition of y. To readily understand the 
physical meaning of the new scaling variable, let us 
consider the asymptotic limit of kminiu, q, E) for a large 
nucleus [i.e., k^^^{y,E) = \y - {E - Emin)]]; it can be 
seen that only when E — Emin does kminiy,E) = \y\, 
in which case B^ = and F"^{y) = f^{y)\ this 



holds only for the deuteron, whereas for a complex 
nucleus E\_^ ^ and E > Emm-, so B^{y) ^ 0, and 
F^{y) 7^ f^{y)- It is therefore the dependence of kmin 
on -E^.j that gives rise to the binding effect [i.e., to 
the relation F^(?/) f^{y)]- This is an unavoidable 
defect of the usual approach to Y-scaling, based on the 
scaling variable y; in fact, the longitudinal momentum 
is very different for weakly bound, shell-model nucleons 
{E\_-^ — 20 MeV) and strongly bound, correlated 
nucleons (-E^_i ~ 50 — 200 MeV), and at large values of 
\y\ the scaling function is not related to the longitudinal 
momentum of strongly bound correlated nucleons, whose 
contributions almost entirely exhaust the behavior of 

13, m, 
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the scaling function. As stressed in Refs. 
to establish a global link between experimental data 
and longitudinal momentum components, one has to 
conceive a scaling variable that could equally well rep- 
resent longitudinal momenta of both weakly bound and 
strongly bound nucleons so that the binding correction 
could be minimized. This can be achieved by adopting a 
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FIG. 1: The average value of E'X-i{k) [Eq. (0] calculated for 
nuclear matter with the spectral function of Ref. (ITl] (open 
dots), and for ^He with the spectral function from the Pisa 
wave functions [3] (full dots). The full lines are obtained with 
the spectral function of the few-nucleon correlation model of 
Ref. 0. 

scaling variable that properly includes the momentum 
dependence of the average excitation energy of {A — 1) 
generated by correlations, namely, 

< E*^_,ik) J P^{k,E\_,)E\_,dE\_, (7) 



where E*, 



E. 



(2) 



E, 



(2) 



- ^ - ^thr^ ^thr - + 2mN - Ma 

being the threshold energy for two-particle emission. We 
have calculated the quantity in Eq. ([7]) using a realistic 
spectral function for nuclear matter and ^He. The results 
are presented in Fig. [1] where they are compared with 
the prediction of the spectral function of the few-nucleon 
correlation (FNC) model of Ref. [l^, according to which 
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where Kc m is the CM momentum of a correlated pair. In 
view of the very good agreement between the FNC model 
and the exact many-body results for nuclear matter and 
"^He, we used the former to calculate < E\_^{k) > for 
nuclei with 3 < A < oo. The values shown in Fig. [T]can 
be interpolated by 



bA - CA|k| 



(9) 



and ca result 
= 37.3 Mel/, 
0.03). Placing 

Ma-i+ < 



where T^v — ( ■\/ mj^ + — tojv ) i and 6^ 
from the CM motion of the pair (bj^M 
CNM = 0.04 and 63 = ~2MMeV, 03 = - 
in Eq. ^ k = yew, ^^^1 = 1 and MX_, 
^A-iik) > — < Egr >, we obtain a fourth-order equa- 
tion for the new scalin g va riable yew, which, in contrast 
to previous work 1^, [isj, [fj], has been solved exactly; 
this, together with the relativistic extension of the defini- 
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FIG. 2: The ratio of the binding correction [Eq. (|6] to 
the scaling function [Eq. ([5])] for ^He (open dots) and ^^C 
(full dots) calculated with the scaling variable y, which does 
not contain any effective excitation energy from SRC (upper 
panel), and with the variable yew, which takes into account 
SRC effects by Eq. @ (lower panel). 

tion of the mean excitation energy, is necessary to extend 
yew to high values. Note that the value of < Egr >, 



fixed by the Koltun sum rule (see Refs. [ll,[ll,[ll), has 
been added to Eq. ^ to counterbalance the effects of 
< E\_-^ > at low yew- For a large nucleus and not too 
large values of yew, one has 



4m% 
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Here, i>a = 
bA+ < E, 
q2 = Ml, + 2vMn - 



Md, 

gr >, q ^ q 

j3 -r ZuMd 



Md = 2mjv - e[1^ - 
- caJ^a and — v^^ — 

For the deuteron yew — 



y = -ql2 + {vd/'2Wd)^/W'^- Amj^ with vd = v + 
Md and invariant mass W|) = vd^ ^ = Md"^ + 
2vMd ~ Q^; for small values of yew, such that 

^^WVCW + "^W - ^^n) + bA- CA\ycw\ « Egr >, 

the variable y, representing the longitudinal momentum 
of a weakly bound nucleon, is recovered. Therefore yew 
effectively takes into account the k dependence of E\_^, 
both at low and high values of yew, and interpolates 
between the correlation and the single-particle regions; 
it can be interpreted as the minimum longitudinal mo- 
mentum of a nucleon that, at high values of yew , has 
removal energy < E'^_^ > and is partner of a correlated 
two-nucleon pair with effective mass Md- 
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FIG. 3: The experimental scaling function (symbols) for *ffe 
^^C, and ^^Fe obtained from the experimental data of Refs 
The upper panel shows F^{q,Y 



y) and the lower 

panel F^{q,Y = yew)- The full, long-dashed, dashed and 
dotted curves represent the longitudinal momentum distribu- 
tions f^{Y) = 27r J^^i n^{k)kdk for '^H, ^He, ^^C and ^^Fe, 
respectively, calculated with realistic wave functions. 



Let us now illustrate the merits of yew-scaling and its 
practical usefulness. The main merit is that, because of 
the definition of yew, binding effects play a minor role, as 
clearly illustrated in Fig. [2l practically kminiq,!^, E) ~ 
\yew\ and B^{q,yew) — 0, with two relevant conse- 
quences: i) to a large extent F^{q,yew) — f^ivcw) [cf. 
Eq. dn])], and ii) as a result of i), one would expect that at 
high values of yew, F^{q, yew) will behave in the same 
way in the deuteron and in complex nuclei, since n^{k) ~ 
CAn^{k) and, accordingly, f^iyew) - CAf°{ycw); at 
low values of yew, in contrast, F^{q,yew) should ex- 
hibit an A dependence generated by the different asymp- 
totic behavior of the nuclear wave functions in configura- 
tion space. This is fully confirmed in Fig. [3]which, more- 
over, also shows that whereas F^{q,y) scales to a quan- 
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FIG. 4: The scaling function F^(Q^ ,ycw) from the lower 
panel of Fig. [3] plotted vs at fixed values of yew He- 
asterisks , ^^C- triangles, Fe- squares). In the right panel the 
data for ^Jfe, ^^C and ^^Fe have been divided by the con- 
stants Ci ~ 2.7, Ci2 = 4.0 and Cse ~ 4.6, respectively. The 
theoretical curves refer to and represent the PWIA results 
(full) and the results that include the FSI (dashed), both ob- 
tained with the AV18 interaction [2Qj. Scaling variables are 
in MeV/c. 



tity that strongly dilfers from the longitudinal momen- 
tum distribution, F^[q, yew) scales exactly to f^[ycw)- 
This is even better demonstrated in Fig. HI where the ef- 
fects of FSI are also illustrated. The left panel shows 
that (i) scaling is violated and approached from the top 
(which is clear signature of the breaking down of the 
PWIA, which has to approach scaling from the bottom 
0) and (ii) the dependence of the scaling violation 
appears to be the same for the deuteron and complex 
nuclei, a fact that has never been demonstrated before 
and represents, in our opinion, a relevant finding. To 
better validate point ii), we have divided F^{Q'^ ,ycw) 
by a constant Ca, such as to obtain F^{Q'^, ycw)/CA — 
(Q^ ,ycw)- The results are shown in the right panel 
of Fig. m it can again be seen that not only at high 
values of Ij/cvkI do all scaling functions scale in A, but, 
more importantly, the constants Ca agree, within the 
statistical errors, with the theoretical predictions of Ref. 
0, as well as with the experimental results on the ra- 
tio R{xBj,Q^) = 2aiixBj,Q^)/Aa^ixBj,Q^) i. The 
main findings of our analysis can be summarized as fol- 
lows: (i) at high values of \ycw\ ( > 200 - 300 MeV/c) 
the scaling function F'^{Q'^,ycw) scales to the one of 
the deuteron, with scaling constants Ca in qualitative 



agreement with theoretical predictions and other types 
of experimental analysis; this kind of A-scaling is en- 
tirely due to the scaling of the momentum distributions, 
n^{k) ~ CAn^{k), at A; > 1.5 — 2 /m^^, which can there- 
fore be investigated by ycwscaling analysis of inclusive 
data, owing to the direct link between the scaling func- 
tion F^{Q'^,ycw) and the longitudinal momentum dis- 
tributions; ii) the FSI has relevant effects on the scal- 
ing functions up to ~ 4 — 5GeV^ but, most impor- 
tantly and surprisingly, it exhibits a similar depen- 
dence in complex nuclei and in the deuteron; this has 
neither been observed nor theoretically predicted previ- 
ously; in a forthcoming paper it will indeed be shown that 
the effects of the FSI on the momentum distribution of 
a correlated nucleon are similar in the deuteron and in a 
complex nucleus (for preliminary results see Ref. p^V 
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